STEADY-STATE FINGERING PATTERNS FOR A PERIODIC 

MUSKAT PROBLEM 



MATS EHRNSTROM, JOACHIM ESCHER, AND BOGDAN-VASILE MATIOC 

Abstract. We study global bifurcation branches consisting of stationary 
solutions of the Muskat problem. It is proved that the steady-state fingering 
patterns blow up as the surface tension increases: we find a threshold value 
for the cell height with the property that below this value the fingers will 
touch the boundaries of the cell when the surface tension approaches a 
finite value from below; otherwise, the maximal slope of the fingers tends to 
infinity. 



1. Introduction 

Proposed in 1934 by Muskat (cf. [11]), the Muskat problem describes the 
evolution of the interface between to immiscible fluids in a porous medium. In 
the recent investigation [6] this problem was studied in a new, periodic, setting 
incorporating gravity, viscosity, and surface tension effects. The current note 
aims at an in-depth description of the stationary solutions found in that work. 

When a heavier viscous fluid rests upon a lighter one, the interface between 
them is in general not stable; depending on the different densities, and the 
surface tension, one expects the upper fluid to, at least partially, sink into 
the lower one, and vice versa. Due to their resemblance to an outstretched 
hand reaching into a viscous fluid, the resulting shapes are often referred to as 
fingering patterns. The investigation of such, in different settings, has brought a 
lot of attention (see, e.g., the pioneering paper [15] and the later investigations 
[4, 10, 13, 14]. 

In [6] smooth branches of stationary, i.e. time-independent, solutions of the 
Muskat problem were found. They are periodic solutions of the Laplace- Young 
equation under a volume constrain (see (2.3)). The Laplace- Young equation is 
also known as the capillarity equation and, subjected to boundary constrains, 
has been studied by many authors (see [7] and the literature therein). 

The solutions we found are all even, but only in a small neighbourhood of the 
trivial solutions can one via linearisation obtain an approximate picture of the 
fingering patterns. This is due to the fact that global bifurcation theorems are 
inherently implicit in nature, and thus have the drawback of not disclosing the 
behaviour of the bifurcation branches away from the bifurcation point. In our 
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FIGURE 1 . The periodic and vertical Hele-Shaw cell. 

present work, we therefore take advantage of the theory for ordinary differential 
equations and certain symmetry properties of the solutions to give a precise 
description of the solutions found in [6]: we show that each global bifurcation 
branch consists entirely of steady-state solutions of minimal period 2ir/l, I £ N, 
and that the symmetric fingers described by the interface i) either approach 
the bottom and the upper boundary of the cell, or ii) display blow-up in the 
C 1 — norm, while the surface tension coefficient tends from below to a finite 
value. 

The plan is as follows. In Section 2 we give the necessary mathematical 
background of the problem, and show that, for stationary solutions, it may be 
reduced to an ordinary differential equation with an additional non-local con- 
straint. The proof of the main result mentioned above is based on the study of 
the odd solutions of this equation, and the one-to-one correspondence between 
the odd and the even solutions thereof. This is done in the Section 3, and there 
we also show that there exist infinitely many global bifurcation branches con- 
sisting of odd solutions of the problem. In addition, we describe the behaviour 
of the steady fingers away from the set of trivial solutions. Finally, it is inter- 
esting to see that the steady-state fingering patterns we obtained correspond to 
certain solutions of the mathematical pendulum. This correspondence is shown 
in the Appendix. 



2. Preliminaries 

Let h > 0, and consider a periodic medium occupying § x [— h,h), with § 
denoting the unit circle. The bottom of this cell is assumed to be impermeable, 
and the pressure on the upper boundary is constantly set to zero. For a function 
/ with ||/|| c( s) < h, let 

r(/(i)) :={(x,f(t,x)) : xGS}, 

be the time-dependent interface separating the wetting phases, and r± := § x 
{i/t} the bottom and the upper boundary of the cell (see Figure 1). We define 
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the fluid domains 



fi-(/(t)) :={(x,y) : -h < y < f(t, x)}, 
«+(/(*)) :={(x,y) : f(t,x)<y<h}, 



and write 



K r(/) 



(1 + fit' 2 

for the signed curvature of the graph The mathematical model can then 

be stated as a two-phase moving-boundary problem, 

' Au± = in n±(f(t)), 



d v u A 



9i 



on r , 



1i+ — U- 



92 on r_, 

7«r(/) + ~ £-)/ 011 r(/(t)), 



(2.1) 



dtf 
/(O) 



on r(/(t)), 



/o 



with t £ [0,T], where we use the subscripts ± to denote the upper and lower 
fluids, respectively. As conventional, g stands for the gravitational constant of 
acceleration, 7 denotes the surface tension at the interface T(f), and g± and 
(j,± are the densities and viscosities of the two fluids, respectively, all of which 
are supposed to be given positive constants. Physically, the potentials u± are 
defined by the relation 

u± ■■=p± + gg±y, 

where p stands for pressure, and y is the height coordinate. Furthermore, the 
functions g\ and §2 are assumed to be known 

9l eC([0,T],h 1+a (S)) and g 2 G C([0,T], h 2+a {S)). 

Given m £ N and a E (0,1), the small Holder spaces h m+a (B>) stand for the 
completions of the class of smooth functions in the Banach spaces C m+a (§). 

The problem consists of finding functions / and u± satisfying (2.1), but it 
can be shown that this may be reduced to a parabolic problem with / as the 
single unknown [6]. Hence, we shall refer to the function / parametrising the 
moving interface between the fluids as a solution of (2.1). 

Well-posedness results. It is shown in [6] that the Muskat problem is, at least 
in a neighbourhood of some flat interface, of parabolic type. This observation 
is true, when considering surface tension effects, independently of the boundary 
data gi and g 2 - On the other hand, when neglecting surface tension, certain 
restrictions must be imposed on the boundary data to ensure parabolicity of 
(2.1). We then have (cf. [6, Theorem 2.1]): 
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Theorem 2.1 (Well-posedness). Let 7 G [0, 00), C\,C2 G K, and assume that 
7>0 or ff (p+-p_)+ci f^- 1 ) <°- ( 2 - 2 ) 

TTien i/iere esisi open neighbourhoods of the zero function Oi C h l+a (E>), i G 
{1,2}, andO C fr 2 + 2si s n (7)+a(S) ; such that for all f G O and 9i G C([0,oo),C; + 
Oj), i = 1,2, i/iere exists T(/o) < T and a unique maximal Holder solution f of 
problem (2.1) on [0,T(/ )) which fulfills f(t) G /or a// 1 G [0,T(/ )). 
7/7 > ; £/ien iwe may choose Oi = /i* +a (§), i G {1, 2}. 

Existence of classical solutions of the Muskat problem, and long-time exis- 
tence for small initial data, can also be found in [8, 16, 17, 18]. The approach in 
[6] yields structural insight into the character of the Muskat and it is suitable for 
studying the stability properties of the steady-state solutions of problem (2.1). 

Steady-state solutions. In the remainder of this paper we assume that g\ = 
and 52 = const, meaning that the mass of both fluids is preserved in time, and 
that the cell contains equal quantities of both fluids. The steady-state solutions 
of (2.1) are then solutions of the problem 

/" f 

7 ,.. , „ 2N o /2 + g{Q+ - g-)f = const, and / / dx = 0. (2.3) 
(1 + fly/* J 8 

Indeed, since / does not depend on time, it follows from uniqueness for the 
Dirichlet-Neumann problem that the potentials n+ and n_ are both constants 
also in the spatial variable, which yields the first equation of (2.3). The second 
relation reflects the earlier mentioned assumption that the cell contains equal 
amounts of both fluids. By induction, we obtain 

Remark 2.2. Any classical solution of (2.3) is smooth. 

We shall refer to the set 

£:={( 7 ,0) : 7 >0} 

as being the trivial branch of solutions of (2.3). Because of the integral con- 
straint in (2.3), the problem (2.3) is in general over-determined. One way to 
approach this difficulty is to determine solution pairs (7, /) G (0, 00) x C 2 (§) of 
(2.3), under the additional, but natural, requirement that ||/||c(S) < ^> meaning 
that the fingers do not touch the lower or upper boundaries of the cell. 

In the situation when the less dense fluid lies on the bottom of the cell, 
i.e. when g + > o_, we find — using the theorem on bifurcation from simple 
eigenvalues due to Crandall and Rabinowitz [3, Theorem 1.7], and the global 
bifurcation theorem due to Rabinowitz [9, Theorem II. 3. 3] — global bifurcation 
branches consisting of even, stationary, finger-shaped solutions of (2.3). More 
precisely, if Cq"^ q (S) denotes the subspace of C 3+Ct (§) consisting of even func- 
tions with integral mean zero, and 

W := {/ G C 3 + Q (S) : \\f\\ m <h}, 

we have (cf. [6, Theorem 6.1 and Theorem 6.3]): 
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Theorem 2.3 (Bifurcation of stationary solutions). Let g± = 0, g<i = const, 
and Qjf. > Q-, 1 < I £ N. The point 

(j l ,0):=(g(g + -g^)/l 2 ,0) (2.4) 

belongs to the closure S of the set of nontrivial solutions of (2.3) in (0, 00) x W. 
Denote by Ci the connected component of S to which (7^,0) belongs. Then Ci 
has, in a neighbourhood 0/(7;, 0), an analytic par ametr is ation (ji, fi) : (—5,5) — > 
(0,oo) x W, 

7^)=7, + ^p^ 2 + 0(, 3 ), 

/,(£) =ecos(/x) + 0(e 2 ), 

as e — > 0. ylny other pair (7, 0), 7 > 0, belongs to a neighbourhood in (0, 00) x W 
on/y trivial solutions of (2.3). 
Furthermore, if e is small and 7 = then fi(e) is an unstable stationary 

solutions of (2.1). 

Theorem 2.3 is obtained by differentiating the first relation of (2.3) and find- 
ing in this way an equation for / only (this is why solutions in Cg^ a (8) are 
considered). It is not difficult to show that if g~ > then (2.3) has only the 
trivial solution / = (see, e.g., [5]). In the paper at hand we show (cf. Re- 
mark 4.3) that this is the case when £>_ < g+ too, as long as the surface tension 
coefficient is large enough. 

3. Odd steady-state fingering solutions 

In this section we consider the odd solutions of (2.3). If / is an odd function 
on §, then / has integral mean and /(0) = f"(0) = 0. Hence, the odd steady 
states of the Muskat problem (2.1) are exactly the odd solutions of the equation 

(1 + /l)3/2 + X f = °> A >°> 

within the set 

W:={/GC 00 (S):||/|| c(§) </ i }. 

Here, the shorthand 

A := 3(e+-Q-) (3 2) 

7 

indicates the character of (3.1) as an eigenvalue problem. Notice that, through- 
out this work, we consider the unstable case when g + > g^, i.e. when the 
heavier fluid occupies the upper part of the membrane. Equation (3.1) admits 
the following scaling property: 

Proposition 3.1. If(^,f) defines a solution of (3.1) through (3.2), then 

(r 2 7 ,r7(0), isn, (3.3) 

is also a solution of (3.1). 
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Proof. Since A is inversely proportional to 7, the result is immediate. □ 

The main result of this work is the following theorem, which states that 
a global bifurcation branch consisting of odd functions of minimal period 2ir 
emanates from the trivial branch of solutions £ at (7^ , 0) , where is defined 
by (2.4). This will later be used to characterise the global bifurcation branches 
of odd solutions which arise at (7^,0),/ > 2, (see Corollary 3.3 below), and in 
Section 4 to describe the global bifurcation branches C[ obtained in Theorem 2.3. 

Recall the definition of the beta function, 

B(x,y):= I f > - x {l-t) v ~ 1 A, Rez,Rey>0. 
Jo 

Theorem 3.2. For each h > 0, there exists 

Aft > A* := ^2 { B (f > l)) 2 > (3-4) 

and corresponding 7^ < 7* defined by (3.2), with the property that the nontrivial 
odd solutions of (3.1) of minimal period 2ir within U coincide with the global 
bifurcation curve 

El := {(7,±/y) : 7 G (7i.7fc)}, 
where the odd function f~ f G C°°(§) is uniquely determined by the parameter 
7 G (7i,7*) if we require that / 7 (0) > 0. Let h* := \j2j A*, and let / 7 denote 
the solution of (3.1) of minimal period 2ir (not necessarily inU). The mapping 
(7i)7*) x§3 (7,1)^ f^( x ) i s smooth, and 
(i) ifh<h*, then 7^ < 7*, and 

\\fl\\c($) = fii^/ 2 ) /~ h as 7 / lh; 
(ii) ifh = h*, then 7^ = 7*, and 

||/ 7 llc(S) = / 7 (tt/2) / h, ||/;|| c(§) = /;(0)/ , oo as 7 /» 7h ; 
(Hi) ifh>h*, then 7^ = 7*, and 

Helices) = / 7 (0) ^ 00 as 7 S lh , 
while sup[- ij7t) ||/ 7 ||c7(S) < h. 

Recall that 7;, 1 < I G N, is the constant defined by (2.4). Combining 
Proposition 3.1 and Theorem 3.2 we conclude: 

Corollary 3.3. Let 2 < I G N. There exists 7^ G (71,7*] with the property 
that 

e, ^{(r^ir 1 /^-)) : 7 G (7l, tu)} 

consists exactly of the nontrivial odd solutions of minimal period 2n/l of (3.1) 
within IA. The alternatives (i) — (Hi) of Theorem 3.2 hold true with the natural 
modifications. The disjoint union 

S 2w : = (U^E,) U (U^ [(7j+i,7 { ) x {0}]) U [(71,00) x {0}] (3.5) 
constitutes all nontrivial 2n— periodic and odd solutions of (3.1) in U. 
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FIGURE 2. A qualitative picture of the 27r-periodic solutions de- 
scribed in Corollary 3.3. 

Remark 3.4. Put differently, Corollary 3.3 states that global bifurcation branches 
consisting of odd solutions emanate from E at 7j, 1 < I € N. Moreover, these 
bifurcation branches are pairwise disjoint. 

Remark 3.5. It is worth mentioning that, for the same 7, we may find 2ir- 
periodic odd solutions of (3.1) of different minimal periods (see Figure 2). Since 
7l,h =7* > 7i) for I large enough, there exists positive integers / G N, such that 

7* 7* 
(Z + l) 2 I 2 

Consequently equation (3.1) possesses a solution which belongs to E; and an- 
other one in £j+i, corresponding to the same 7. 



In order to prove Theorem 3.2 we need some preliminary results. 
Proposition 3.6. Let A > and a £R be given. The initial-value problem 

f" 

(1 + ^2)3/2 = 071 § ' 

/(0) = 0, ( 3 - 6 ) 
f(0) = a 

possesses a unique classical solution f\ a . The solution is odd and periodic in 
x, and smooth as a map 

(0,oo) x M x fx, a (x). 
Proof. Setting g := /', we rewrite (3.6) as an initial value problem for the pair 

where F : R 2 — > R 2 is defined by 

F ( = ( "A/(l + 5 2 ) 3/2 

Since F is smooth, there is a unique and smooth solution of (3.7), defined on 
a maximal interval [0,T); if T < 00, then the solution blows up, meaning that 
sup[ ,T) \(f,g)\ = 00 (cf. [1]). 
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Notice that if / is an odd solution of (3.1) with slope /'(0) = a > 0, then 
— / is also an odd solution of (3.1) with slope —a. Without loss of generality 
we may therefore restrict our attention to solutions of (3.1) with nonnegative 
slope at x = 0. Clearly, the solution of (3.6) with slope a = is / = 0. 

Suppose now that a > 0. We prove that there exists a positive constant 9 x ^ a 
such that /' > on [0,0 X a ) and f'(0x,a) = 0. Indeed, assuming the contrary, 
we obtain in view of /'(0) = a > 0, that /' > on [0, T). 

On the one hand, if T = oo, we infer from (3.1) that 



= /"(x) + A/(s)(l + (/'(s)) , 
Integration yields that 

/'(x)</'(l)-A/(l)(x 



> /"(x)+A/(l), for all x>l. 



1) 



-oo 



as 



oo, 



which contradicts our assumption. 

On the other hand, if T < oo, then either supp T ^ f = oo or supj ^ f = oo, 
the latter case being excluded by the fact that /' is decreasing for positive /. 
If sup[ 0jT ) / = oo, we multiply (3.1) by — /' and integrate over [0, x] to obtain 
that 

1 1 , Xf(x) 



< x < T. 



(3. 



(l + / /2 (x))V2 (l + a 2)l/2 2 

Letting x — > T, we obtain the desired contradiction. Consequently, there exists 
a unique 6\^ a > 0, such that /' > on [0,6\ tO[ ) and f'(0x,a) — 0. 

It can be easily seen that / extends to an odd function of minimal period 
T\ )Ce := ^0\ a . Indeed, we see that 

( 



fx,, 



m, 

f(26 Xta -x), 
-/(s-20 A , a ), 
-f(4Bx, a -x), 



< x < A , a , 

0\,a <X< 28 X , a , 

26 Ka <x< 36 XtCt , 

3^A,a < X < 40 A , a , 



has an odd and T A q ,— periodic extension on the whole of ' 



(3.9) 



□ 



We now explicitly determine the minimal period, called T Xa , of the solution 
fx, a °f (3.6). In order to simplify calculations, we put 



1 



Q" 



From relation (3.8), we find for x = 6 x>a that the maximum of f Xa is 



fx,a(0X,a) = V2A-M1-/3)- 

We also infer from the same relation that 



(3.10) 



f'xM 



+ 



1 for all x G [0, 6. 



7,aJ 
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Dividing this equality by its right-hand side, we find that 

'tJU+^MY-!] dx , 



'X.i 







and the variable substitution f(x) = s yields 

ff\,a(Q\,a) 

6\,a = / 
J 




1 ds. 



Finally, setting r := s/ f\ a (9\ a ), we obtain in virtue of (3.10) that 

- Ar . v-er+p ( 3.n) 

for all a, A > 0. Since a i— )• /3 is smooth, we may extend 9\ a continuously to 
(0,oo) x [0, oo). More precisely, we state: 

Lemma 3.7. The function 9\ a defined in (3.11), 

(0,oo) x [0,oo) 9 (A, a) (->■ 9 x , a G (0,oo), 

is smooth, and strictly decreasing with respect to both A and a. Moreover 1 

9 Xt0 = and lira A)Q = B f £Y (3.12) 

2VA 2V2A V 4 2/ 

Proof. The integral on the right-hand side of (3.11) exists because the singu- 
larity behaves like (1 — r) -1 / 2 as r — > 1. Therefrom, the regularity assertion is 
clear. Let us now show that 9\ a is strictly decreasing with respect to a. To 
this aim we fix r G (0, 1) and define the function 

g T (a) := (l-/3)r 2 + /3 

= (1 - (1 + a 2 y l ' 2 )T 2 + (1 + a 2 )" 1 / 2 , a > 0. 

Since 

^ f 1 1 <?r(a) 



dr, 



A7 Vl-r 2 71 + ^(0) 
we see that 9\ a is strictly decreasing with respect to A, and it suffices to show 
that the mapping [[0, 00) 3a4 g T (a)(l + g T (a)) -1 / 2 ] has a negative derivative 
for all r G (0, 1), a > 0. Indeed, since for such a and r we have that 

dg T o ®P a 

-W = 1 ~ T >0 ^ = -(l + a 2 )3/ 2 <0 ' 



^Let 7 > be fixed. Recall that, given a > 0, the value Ta, q = 4#A,a denotes the minimal 
period of the solution f\ tCt of (3.6) and that the latter problem possesses the trivial solution 
fx,o = if a = 0. Having said this, it is clear that the value #a,o = 7r/(2\/A) is not related to 
the trivial solution /a,o = 0, but is just the limit of 6x, a as a \ 0. 
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it follows from the chain rule that 

d ( 9t(p) \ 9t(«)( 2 + 5r(")) 



da\l + g T {a)) 2(1 + g T {a)fl 2 



< 0, for all r G (0, 1) 



In view of that g T (0) = 1, the first equality in (3.12) follows. Taking into 
consideration that lim Q ,_ >00 g T (a) = r 2 , we obtain that 

lim 9 X a = */? , T " = — ^= s-y\l - sy 1 / 2 ds 

a^oo V A J 2V2X Jo 

1 ..b(U 



2V2X V4'2, 

This completes the proof. □ 

Recall that we are interested in determining the solutions of (3.1) which 
are not only odd, but also of minimal period 2tt. Thus, we are interested in 
determining the set of A and a such that 9\ : a = ir/2. The following lemma 
provides an answer in terms of a function A 1— > a. 

Lemma 3.8. Let A* be the constant defined by the relation (3.4). If X $ (A*, 1] 
then 9\ tOC 7^ vr/2, but given A G (A*, 1] there exists a unique a(A) G [0, 00) such 
that #A,a(A) = 7r /2- The mapping 

(A*,l] 9 A H> a(A) G [0,oo) 

is smooth, bijective, and decreasing. 

Proof. In view of Lemma 3.7, we have that 9\^ a = ir/2 if and only if 

1 / 3 1 \ ^ 77 

which is equivalent to that A G (A*,l]. Since [(A, a) 1— > 0\ tOl ] is smooth and 
d a 9\,a < for a > 0, we infer from the implicit function theorem that [A 1— > 
a(A)] is smooth as well. Then 

= ^y0\ a (\) = d\6 x , a + {d a 0x, a ) a' (A), 

so that c/(A) < for A G (A*, 1) in view of Lemma 3.7. From (3.12) we infer 
that a(l) = and lim^A, Q (A) = 00. □ 

With these preparations done, we come to the proof of the main result as 
stated in Theorem 3.2: 

Proof of Theorem 3.2. It follows from Proposition 3.6 and Lemma 3.8 that the 
odd solutions of (3.1) of minimal period 2tt coincide with the set 

{(A,±/ A ) : A G (A*, 1]}, 

where we simply write fx '■= /a,o(A)- m order for those solutions to be physically 
realistic, we still have to require that f\ G U, i.e. that ||/||c(S) < h. The 
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FIGURE 3. Steady-states fx on Si, A* < A2 < Ai < 1. 



maximum of \f\\ is achieved at x = #a,«(A) = 7r /2, and we infer from (3.10) that 
A £ (A*, 1] must additionally satisfy 

/a(tt/2) = yJvFi (l-(l + a 2 (A))^V2) < h . (3.13) 

Let us first assume that h 2 < 2A" 1 . Since, in view of Lemma 3.8, 

2A- 1 (l-(l + a 2 (A))~ 1 / 2 ) / A _ A . 2A; 1 

we find a unique A^ > A* with the property that 

max I f Xh I = ^2A h 1 (l-(l + a2(A,,))-i/2) = h . 

By recalling that A = g(g+ — Q-)/"f, we infer the main part of the theorem, 
including (i). 

If instead h 2 = 2A" 1 , then fx € U for all A G (A*, 1] and Lemma 3.8 implies 
that 

\h\ /V+A, \/VK, while f x (0) = a(A) 00. 
We have thus shown that (ii) is valid, and assertion (Hi) follows similarly. □ 

Remark 3.9. Differentiating relation (3.8) with respect to A and using a max- 
imum principle argument, shows that f\ ± > f\ 2 on (0, ir) provided A* < Ai < 
A2 < 1. The evolution of the solution fx with respect to A € (A*, 1] is pictured 
in Figure 3. 



Proof of Corollary 3.3. The constant 7^ is defined similarly to 7^, and ensures 
that Z -1 / 7 (Z-) remain in IA for all 7 G [litJih) ( see the proof of Theorem 3.2). 
If I is large enough relation (3.10) shows that 7^ = 7*. The relation (3.5) 
now follows in view of that the global bifurcation branch consists exactly of 
functions with minimal period 2n/l. □ 
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4. Description of the bifurcation branches Q 
Let us return to the setting of Theorem 2.3. Define 

% ■= {(7, fy(- + tt/20) : (7, h) £ U {(71, 0)}, 1 < I £ N. 

Since the functions / 7 are odd (cf. (3.9)), the smooth curve £/ consists of even 
functions. Hence, it must be a subset of the maximal connected component 
of S to which (7j,0) belongs, i.e. Y>i C C\. We shall prove that the converse 
is also true. This means that the branches C% consist, with the exception of 
the trivial solution (77,0), exactly of even functions with minimal period 2n/l. 
Theorem 3.2 may be then used to describe the global bifurcation branches Ci. 

Theorem 4.1. Given I £ N, I > 1, we have that Q = £j. 

Since the even solutions of (2.3) near (7^,0) lay either on the trivial curve 
E or on C\, we conclude that and C\ coincide in a small neighbourhood 
of (7;,0). Hence, at least in small neighbourhood of (7^,0), the (seemingly) 
arbitrary constant in (2.3) is zero. Even more holds: 

Lemma 4.2. Let (7,/) £ (0, 00) x W be a solution of (2.3). We then have 

Moreover, if2n/l,l > 1, is the minimal period of f, then /(• — tt/21) is an odd 
solution of (2.3). 

Proof. Assume by contradiction that we would find a solution (7, /) £ (0, 00) x 
W of (2.3) such that 

f" 

(1 + //2)3/2 + X f = C 111 S > 

with a constant c ^ 0. Let g := / — c/A. Then g is an even solution of (3.1), but 
g has no longer integral mean equal to 0. Since </(0) = 0, it must hold g(0) 7^ 0. 
Otherwise, g = 0, meaning that / = c/A, which contradicts / £ VV. 

We may, without loss of generality, assume that g(0) > 0. There exists a 
positive time T c > such that g > on [0, T c ) and g(T c ) = 0. Indeed, if this is 
not the case, we infer from (3.1) that g" < 0. This is in contradiction with the 
fact that g is periodic. The function g is a 4T C — periodic function on M, since it 
must hold that 

g(x), < x < T c , 
-g(2T c -x), T c <x<2T c , 



g(x) 



(4.1) 

-g(x-2T c ), 2T c <x<3T c , 



g(4T e -x), 3T c <x<4T c . 

Moreover, g is 2tt— periodic, so that T c = ir/(2k), for some k £ N. We conclude 
that g has integral mean zero, which is in contradiction with / £ W and c 7^ 0. 
Thus, / must solve (3.1). The relation (4.1) then holds also for /, provided that 
T c = n/ 21. This completes the proof. □ 
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In virtue of Corollary 3.3. Proposition 3.6, and the proof of Lemma 4.2 we 
conclude: 

Remark 4.3. The solutions of (2.3) are, up to translation, odd. Moreover, 
problem (2.3) has no solutions (7, /) G (0, 00) x U if 7 > 7/,. 

With this preparation done, the proof of Theorem 4.1 is immediate. 

Proof of Theorem 4-1- Lemma 4.2 shows that the mapping 

U£i (£, U {(7, , 0)}) -+S, X, U {(7,, 0)} 9 (7, /) -> (7, /(• + tt/20), 

is one-to-one and onto. Recall that S is the closure of the set of nontrivial 
solutions of (2.3) in (0, 00) x W, and the union (£j U {(7;, 0)}) is, in virtue 
of Corollary 3.3, disjoint. The image of £/ U {(t~j, 0)} under this mapping is 
hence C\ = S;. □ 

Appendix 

It is clear form Proposition 3.1 and Corrolary 3.3 that the T— periodic solu- 
tions of (2.3) are exactly the elements of the set 

St =={/(!•) =/^4 

The volume constrain in (2.3) for T— periodic solutions reads as f^fdx = 
0. In view of Remark 4.3 we show that each T— periodic solution / of (2.3) 
corresponds to a unique function 9 describing the evolution of a mathematical 
pendulum (or the bending of an elastic rod): 

#" + Asin(#) = 0, (4.2) 

cf. also [12]. We set A := g/l, with I denoting the length of the pendulum. 
We refer to [1, 2] for the deduction of (4.2). Herein, of interest are only the 
solutions of (4.2) which satisfy 

|0| < vr/2 and 9(0) = 0. (4.3) 

Particularly, solutions of (4.2)-(4.3) are odd. We now state: 

Theorem 4.4. There exists a one-to-one correspondence between the even so- 
lutions f of (2.3) and the odd solutions 9 of (4.2) -(4.3). 

Given s G R, 9(s) is the angle between the tangent to T(f) at z(s) and the 
Ox-axis, with z : R — > M a parametrisation ofT(f) by the arc length. 

Proof. Take first / to be an even solution of (2.3). We define the function 
p : R -)• R by 



p(x) 



Jo 



This mapping is bijective, and let z denote its inverse. Let 9 : R — > R be given 
by 

9(s) = arctan/'(z(s)), s G R. (4.4) 
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If L = p(T), then 9 is L— periodic. Indeed, we have that 

rz(s+L) rz(s)+T 

L = p(z(s + L)) - p(z(s)) = ^/TTTHf)dt= Vl + f' 2 (t) dt, 

J z(s) J z(s) 

hence z{s + L) = z(s) + T. Being a composition of odd functions, 9 is also odd. 
Given s£l, 

6 ' {S) = 1 + /*(*(l)/ (a) = ~ X f^ z '^ Vl + f' 2 (*(s)) = -*/(*(«)), 
due to 

= ^)) = (l + /'2(z( S )))l/2- (45) 



Hence 



?"( S ) = A/'CsOOyOO = -A 7=£^= = = -Asin(0( S )), 



and 9 is a L— periodic solution of (4.2)-(4.3). 

Conversely, given solution of (4.2)-(4.3), we let z : R — >• R be the function 
defined by 

z(s) := / cos{9(t))dt, s£l, 



and write p for its inverse. Setting T := J cos(0(i)) (it, we have that p(x + T) 
L + for all x € R. It can be then easily verified that / : R — > R, 

: = + / tan(0(p(t))) (it, iGR, 



A 

is a T— periodic and even solution of (2.3). □ 

With this observation, our result stated in Theorem 3.2 rewrites for the math- 
ematical pendulum equation as follows: 

Corollary 4.5. There exists a smooth curve 9\, A G (A*,l], consisting of 
L\— periodic solutions of (4.2)-(4.3) with the property that 

sup|#| = arctan(a(A))) /*\-+\,, 7r /2. 

Remark 4.6. It is worth noticing that the period Lx of these solutions is strictly 
decreasing with respect to A. Indeed, it holds that 

2 r' 2 d9 

A y/X J-tt/2 yj\ - sin 2 (arctan(a(A)) /2) sin 2 (6>) ' 

and since a decreases with respect to A we obtain the desired conclusion. 
Though L\ can be calculated in terms of elliptic integrals, it is in general dif- 
ficult to specify for which solution 9 of (4.2)-(4.3) of period L > it holds 
that z(L) = 2ir, that is the corresponding solution of (2.3) has period 2tt. Fur- 
thermore, the result stated in Corollary 4.5 can not be obtained via standard 
bifurcation theorems, since the period of Lx must decrease with respect to A. 
These facts serve as a motivation for our approach. 
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